Theory of periodic swarming of bacteria: application to Proteus mirabilis 
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The periodic swarming of bacteria is one of the simplest examples for pattern formation produced 
by the self-organized collective behavior of a large number of organisms. In the spectacular colonies 
of Proteus mirabilis (the most common species exhibiting this type of growth) a series of concentric 
rings are developed as the bacteria multiply and swarm following a scenario periodically repeating 
itself. We have developed a theoretical description for this process in order to get a deeper insight 
into some of the typical processes governing the phenomena in systems of many interacting living 
units. 

Our approach is based on simple assumptions directly related to the latest experimental obser- 
vations on colony formation under various conditions. The corresponding one-dimensional model 
consists of two coupled differential equations investigated here both by numerical integrations and 
by analysing the various expressions obtained from these equations using a few natural assump- 
tions about the parameters of the model. We have determined the phase diagram corresponding 
to systems exhibiting periodic swarming and discuss in detail how the various stages of the colony 
development can be interpreted in our framework. We point out that all of our theoretical results 
are in excellent agreement with the complete set of available observations. Thus, the present study 
represents one of the few examples, where self-organized biological pattern formation is understood 
within a relatively simpe theoretical approach leading to results and predictions fully compatible 
with experiments. 



I. INTRODUCTION 



To gain an insight into the development and dynamics 
of various multicellular assemblies, we must understand 
how cellular interactions build up the structure and re- 
sult in certain functions at the macroscopic, multicellular 
level. Microorganism colonies are one of the simplest sys- 
tems consisting of many interacting cells and exhibiting 
a non-trivial macroscopic behavior. Therefore, a num- 
ber of recent studies have focussed on experimental and 
theoretical aspects of colony formation and the related 
collective behavior of microorganisms fljj^]. 

The swarming cycles exhibited by many bacterial 
species, notably Proteus (P.) mirabilis, have been known 
for over a century || . When Proteus cells are inoculated 
on the surface of a suitable hard agar medium, they grow 
as short "vegetative" rods. After a certain time, however, 
cells start to differentiate at the colony margin into long 
"swarmer" cells possessing up to 50 times more flagella 
per unit cell surface area. These swarmer cells migrate 
rapidly away from the colony until they stop and revert 
by a series of cell fissions into the vegetative cell form, 
in a process termed consolidation. The resulting vegeta- 
tive cells grow normally for a time then swarmer cell dif- 
ferentiation is initiated in the outermost zone (terrace), 
and the process continues in periodic cycles resulting in a 
colony with concentric zonation depicted in Fig ^. Sim- 



ilar cyclic behavior has been observed in an increasing 
number of Gram-negative and Gram-positive genera in- 
cluding Proteus, Vibrio, Serratia, Bacillus and Clostrid- 
ium (for reviews see 

The reproducibility and regularity of swarming cycles 
together with the finding that its occurrence is not lim- 
ited to a single species suggest that periodic swarm- 
ing phenomena can be understood and quantitatively 
explained on the basis of mathematical models. In 
this manuscript we first give an overview of the rele- 
vant experimental findings related to the swarming of P. 
mirabilis, then construct a simple model with two limit 
densities. We then investigate the behavior of the model 
as a function of the control parameters and compare it 
to experimental results. 



II. OVERVIEW OF EXPERIMENTAL FINDINGS 

a. Differentiation. As reviewed in |^||, the differen- 
tiation of vegetative cells is accompanied by specific bio- 
chemical changes. Swarmer cells enhance the synthesis 
of flagellar proteins, extracellular polysaccharides, pro- 
teases and virulence factors, while exhibit reduced over- 
all protein and nucleic acid synthesis and oxygen uptake. 
These findings may be explained by arguing that the pro- 
duction and operation of flagella is expensive and may 
require the repression of non-essential biosynthetic path- 
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ways. The largely (10-30 fold) elongated swarmer cells 
develop by a specific inhibition of cell fission which seems 
not affecting the doubling time of the cell mass or DNA. 




FIG. 1. Typical Proteus mirabilis colony. It was grown on 
the surface of a 2.0% agar substrate for two days. The inner 
diameter of petri-dish is 8.8 cm. Gray shades are proportional 
to the cell density: the cyclic modulation is apparent. 

The differentiation process is initiated by a number 
of external stimuli including specific signalling molecules 
and physico-chemical parameters of the environment. As 
an example for the latter, the viscosity of the surrounding 
medium is presumably sensed by the hampered rotation 
of the flagella Neither the signal molecules that 

initiate the differentiation nor the involved intracellular 
signalling pathways are identified yet, but a correspond- 
ing transmembrane receptor has been found recently Q . 
The structure of this receptor, together with other find- 
ings reviewed i n jpfl , suggest a 1 quorum sensing' regula- 
tory pathway |10| characteristic for many, cell density 
dependent collective bacterial behaviors like sporulation, 
luminescence, production of antibiotics or virulence fac- 
tors Q. 

b. Migration of swarmer cells. It is well established 
|^| that swarmer cell migration does not require exoge- 
nous nutrient sources, since swarmer cells replated onto 
media devoid of nutrients continue to migrate normally. 
The ability of migration depends on the local swarmer 
cell density as isolated single swarmer cells cannot move, 
while a group of them can. It was also demonstrated |l2j ] 
that the mechanism by which bacteria swarm outwards 
involved neither repulsive nor attractive chemotaxis. The 
typical swimming velocity (i.e., in liquid environment) 
of swarmer cells is rj lOOmm/h Jl3| which is also their 
maximal swarming speed and rate of colony expansion on 
soft agar plates H . In the usual experimental conditions 
for investigating swarming colony formation the front ad- 
vances with a speed of 0.5 — 10 mm/h H,^ [l5ljl3 ] . Unfor- 



tunately, in these cases there is no information available 
on the velocity of individual swarmer cells, but it must be 
between the colony expansion speed and the swimming 
velocity. 

c. Consolidation. The molecular mechanisms of con- 
solidation, i.e., the downregulation of the gene activity 
responsible for swarming behavior pi is even less known 
than that of differentiation. If the swarming motility uti- 
lizes intracellular energy reserves as has been suggested 
H , then swarmer cells must have a finite lifetime. In ad- 
dition to the septation of swarmer cells taking place at 
the outermost terrace, inside the colony the differentia- 
tion process, i.e., the supply of fresh swarmer cells must 
also be shut off. The cessation of swarmer cell production 
does not seem to be due to severe nutrient depletion since 
vegetative cells keep growing (although with decreasing 
growth rate) well inside the colony for many hours after 
the last swarmer cells were produced in that region Q . 

d. Colony formation. The cycle time (total length of 
the migration and consolidation periods) has been found 
|p^| to be rather stable (« Ah) for a wide range of nu- 
trient and agar content of the medium. The size of the 
terraces and the duration of the migration phases were 
strongly influenced (up to an order of magnitude, and 
up to a factor of 3, respectively) by the agar hardness. 
The nutrient concentration did not have an observable 
effect on these quantities from 0.01% up to 1% glucose 
concentrations. There was, however a remarkable posi- 
tive correlation between the cycle time and the doubling 
time (ranging from 0.7h up to 1.8h) of the cells. 

Two interfacing colonies inoculated with a time dif- 
ference of a few hours and therefore being in different 
phases of the migration-consolidation cycle, were found 
to maintain their characteristic phases Jl4||l5| ■ Thus, the 
control of the swarming cycle must be sufficiently local. 

e. Cycle rescheduling. A few experiments investi- 
gated the cell density dependence of the duration of qui- 
escent growth [lag phase) prior to the first migration 
phase. These studies clearly revealed that vegetative 
cells have to reach a threshold density to initiate swarm- 
ing fl4| , |l5|| , in accord with the suggested quorum sensing 
molecular pathway of the initiation of swarmer cell dif- 
ferentiation. 

Agar cutting experiments demonstrated that a cut in- 
side the inner terraces does not influence the swarming 
activity |]l5| . However, when the cut has been made just 
behind the swarming front, the duration of the swarming 
phase was shortened and consolidation was lengthened 
by up to 40% jn|. 

Even more interestingly, mechanical mixing of the cell 
populations before the expected beginning of consolida- 
tion expands the duration of the swarming phase consid- 
erably, by up to 50% p3{ . This finding, together with 
replica-printing experiments Jl3[ | demonstrates that at 
the beginning of the consolidation phase still a large pool 
of swarmer cells exists and seems to be "trapped" at the 
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rear of the outermost terrace. 



III. THE MODEL 

Taking into account the above described experimental 
findings, here we construct a model which is capable of 
explaining most of the observed features of colony expan- 
sion through swarming cycles. 

(i) The model is based on the vegetative and swarmer cell 
population densities only, denoted by p° and p*, respec- 
tively. These values are defined on the basis of cell mass 
instead of cell number, therefore one unit of swarmer cells 
is transformed into one unit of vegetative cells during 
consolidation. 

(ii) Vegetative cells grow and divide with a constant rate 
ro « lh _1 ]T^ , p^ |. This will later allow us to establish a 
direct correspondence between cell density increase and 
elapsed time. 

(iii) Usually, swarmer cell differentiation is initiated when 
the local density of the vegetative cells exceeds a thresh- 
old value (p° mm w 1CT 2 cells/^m 2 [0Jl|). (Prior the 
first swarming phase experiments indicate the presence of 
an extra time period tt which is probably associated with 
the biochemical changes required to develop the ability 
of the swarming transition. This effect is present only at 
the seeding of the colony, thus tt = otherwise.) When 
p° = p° min at time to, some of the vegetative cells en- 
ter the differentiation process, modeled by introducing 
a rate r. Since the biomass production rate is assumed 
to be unchanged during the differentiation process, the 
rate of producing new vegetative cells is ro — r and the 
differentiating cells elongate with the normal growth rate 
ro- 

(iv) The full development of swarmer cells, i.e., a typical 
20- fold increase in length needs a time (td ~ ln20/ro ~ 
3h) comparable with, or even longer than the length of a 
consolidation period, hence can not be neglected. There- 
fore, the first "real" swarmer cells, which are able to move 
appear only at to + t<i ■ 

(v) The production of swarmer cells is limited in time 
and the length t of the time interval during which vege- 
tative cells can enter the differentiation process is another 
phenomenological parameter of our model. 

As we here focus on the periodicity of colony expan- 
sion, we do not consider what happens in the densely 
populated regions after to + r. Specifically, in our model 
any activity of the vegetative cells ceases in these parts 
of the system. 

(vi) Swarmer cells can migrate only if their density ex- 
ceeds a threshold density p* min - Above that threshold, 
swarmer cells are assumed to move randomly with a dif- 
fusion constant Dq. The finite lifetime of swarmer cells 
is incorporated into the model through a constant rate 
(r*) decay. 



Unfortunately there are no good estimates on these 
parameters in the literature. According to recent exper- 
imental observations ]l7j , p* min is less than 10% of the 
value of p° prior the beginning of the swarming phase, 
i.e., 0.1p^ nin e r ° td w 10~ 2 cells/^im 2 . D Q may be esti- 
mated as v^t p /2 with vq being the typical speed of swarm- 
ing cells and t p being the persistence time of their motion. 
As vq is approximately 30mm/h (see Sec II. b) and t p is in 
the order of minutes, Dq is estimated to be in the order 
of 10 mm 2 /h. 

The above considerations lead to the following set of 
equations 

p°(t) = r oP °(t) + r(t)~r*(t) 

p *(t) = -r(i) +r*(t - t d )e rotd + VD{p*)v P * (i) 

where the consolidation (r) and differentiation (r*) 
terms are given by 

T = r*p* and T* =r(p°,t)p°, (2) 



(a) (r) 




FIG. 2. Schematic representation of the model for swarm- 
ing colony formation, (a) The two basic quantities are the veg- 
etative and swarmer cell densities, which can be transformed 
into each other, and changed by proliferation and motility. 
The dotted lines represent regulatory (threshold) effects: The 
rate of differentiation is assumed to be dependent on the veg- 
etative cell density, and the motility of swarmer cells is also 
determined by their local density, (b) Notations and typical 
time courses of vegetative and swarmer cell densities prior 
the beginning of the first migration phase. If t < to the den- 
sity p (or p°) grows with a constant rate ro. After reaching 
the density threshold Pmim vegetative cells keep growing only 
with a rate ro — r, but the total density (p) of the vegetative 
cells and the differentiating swarmer cells still grows with a 
rate ro. After the time required for the full elongation of a 
swarmer cell (when p reaches p m %n), p* becomes positive, and 
for r* = asymptotically would grow with a rate r. 
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Eqs. ([j]) can be significantly further simplified by mak- 
ing use of the possibility to measure elapsed time with 
the increase in p° . In particular, neglecting correction 
terms related to consolidation in areas where p° > p min , 
i.e., where primarily differentiation takes place, we can 
cast r*(t - t d )e rotd in the form 



rotd _ 



rp°(t -t d )e 



' Pmin 

e 



e (r -r)(t-t d -t ) e r t d 



fo(t-to) p -r(t-to-td) — rp°(t)e rtd 



(3) 



for t > to + td — t\. Let us introduce a transformed 
population density p as 



P(t) 



P°(t) 
p°{t)e rtd 



for 



p < p - 

r rmin 

for p°>pg 
for p > po 



i.e., for t < to 
and t < ti 
and t > ti 



(4) 



with a = 10 providing a rather steep, but continuous 
crossover. 
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which is in fact the total density of the vegetative and the 
differentiating, but not yet fully differentiated swarmer 
cells (see Fig. ||). With this notation, using (^) and sim- 
ilar considerations for to < t < to + td, Eqs. (0) can be 
written into a simple, not retarded form 
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p = r Q p + r*p* -r(p)p 
p* = -r*p* + r(p)p + VD{p*)Vp\ 



where 



r(p) 



r for p m i n < p < p rnax 
otherwise 



(5) 



(6) 



with p m 

_(r -r)r 



f» 10 1 cells/^m 2 and p ma x{r) = 
If the T <^ (r — r)p condition does not 



p rtd n 

e P' 



hold in the [to, *i] time interval, then p m i n must be also 
treated as a dynamical variable. This case will not be 
considered here. In the following we use p for the char- 
acterization of vegetative cell density. 



IV. RESULTS 



A. Numerical method in ID 
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FIG. 3. Time development of the model obtained by nu- 
merical integration of the equations, starting from a local- 
ized "inoculum" at t = 0, x = 0. The continuous line rep- 
resent the colony boundary (maximal value of x for which 
p(x) + p*(x) > 0). The filled gray and black areas are regions 
where swarmer cells are motile, and where swarmer cells are 
produced, respectively. For p ma x = 1.3, pmin = 0.01, r — 0.3 
and r* = 1.0 the expansion of the system is clearly periodic 
(a). If we increase the production of swarmer cells by increas- 
ing pmax to 2.0 then the periodicity is gradually lost and a 
continuous expansion takes place (b). 



The model defined through Eqs. (||) has the following 
seven parameters: the rates ro, r, r* , threshold densities 
Pmin, Pmax (or t), p* min and the diffusivity D a . How- 
ever, this number can be reduced to four by casting the 
equations in a dimensionless form using 1/ro ~ lh as 
time unit, p m i n ~ 0.1 cells/^m 2 as density unit and x$ = 

Dq j To ~ 3mm as the unit length. The resulting control 
parameters are r/r , r*/r a , p ma x/ ' Pmin = exp[(r - r)r] 
and p m i n / Pmin- To obtain continuous density profiles, 
the step-function dependence of D on p* was replaced by 



D(p*) = 



D 



1 + tanh 2a- 



(J) 



Representative examples for the time development of 
the model are shown in Fig. ^. The production of 
swarmer cells is localized, and determined by the density 
profile of vegetative cells at the end of migration peri- 
ods. In this particular model p(x) is decreasing towards 
the colony edge, therefore in the migration phases the 
source of swarmer cells is moving outwards. The front 
of swarmer cells is expanding from the inside of the last 
terrace. Because of the decay term T, cells become non- 
motile first at the colony edge. 
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B. Phase diagram 

Each of the dimensionless control parameters can have 
an important effect on the dynamics of the system. As an 
example, if the duration r of swarmer cell production is 
increased, then the consecutive swarming cycles are not 
separated and a continuous expansion takes place with 
damped oscillations (Fig. ||b). To map the behavior of 
the system as a function of the control parameters, the 
following procedure was applied. Migration periods were 
identifed by requiring max x p*(x) > p* min - For a given 
set of parameters we determined the lengths {ti} of the 
consecutive migration periods, and the system was clas- 
sified as periodic if the three largest values of the set 
{ti} were the same within 20%. Otherwise, the expan- 
sion was classified as continuous as long as max^jti} was 
large enough: comparable with the total duration of the 
simulated expansion. 
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FIG. 4. Phase diagram of the system as the function of 
cumulative swarmer cell production density P and migration 
threshold p* min for various values of swarmer cell decay rate 
r* . If the production is high or the motility threshold is low 
enough then continuous expansion can be observed. On the 
other hand, if the production is too low or the motility thresh- 
old is too high then no expansion takes place. In an intermedi- 
ate regime periodic growth can be observed. The boundaries 
of this parameter regime are plotted for r* /ro = 1.0 (thick 
continuous line), 0.4 (thick dashed line) and 0.1 (thick dot- 
ted line). The thin continuous line represent an approximate 
upper bound (^) for cyclic colony expansion. The insert 
demonstrates that the fourth parameter of the model, r, is 
irrelevant: for r* /ro = 1 and r/ro = 0.01 (□), 0.3 (o), 0.5 
(A) and various values of p m ax and p* m i n the type of colony 
expansion was classified. Open symbols correspond to cyclic 
growth, filled symbols to continuous growth while dots denote 
no expansion. Note that the corresponding regions completely 
overlap irrespectively of the value of r. 

The behavior of the model is summarized in Fig. ^, 
where the boundaries of the various regimes are plotted 
for three different values of r* /r . We found, that r and 



Pmax can be combined into one relevant parameter, the 
swarmer cell production density, as 



P = 



T*{x,t)dt = 



ro 



-(pr, 



(8) 



which quantity does not depend on the choice of position 
x. As the insert demonstrates, for a given P, the actual 
values of r or p max are irrelevant to this kind of classifi- 
cation in the parameter regime investigated. The general 
structure of the phase diagram was found to be similar for 
various values of r* . For large enough P or low enough 
Pmin a continuous expansion takes place, while for too 
small P or large p^m the expansion of the system is 
finite. For intermediate values of these parameters an 
oscillating growth develops exhibiting well distinguish- 
able consolidation and migration phases. As the lifetime 
of the swarmer cells is increased, the parameter regime, 
in which periodic behavior is exhibited, is shrinked and 
moved towards lower P values. 

One can easily estimate the position of the boundary 
of the non growing phase based on that (i) the width w 
of the terraces is small (this assumption is justified later, 
in Fig. ^.), thus (ii) the time required for the diffusive ex- 
pansion of the swarmer cells is much shorter than their 
lifetime, which, in turn, is (iii) shorter than the duration 
of a swarming cycle: r* /tq ~ 1 . The amount of swarmer 
cells produced in one period is Pw. Neglecting the decay 
during expansion, the width w' of the next, new terrace 
can be determined from the conservation of cell number 
as 



2p: 



,w 



w(P + p* - p* min ), 



(9) 



where p* denotes the swarmer cell density remaining from 
the previous swarming cycle and the symmetric expan- 
sion of the released swarmers was also taken into account. 
To achieve a sustainable growth w' > w is required, re- 
sulting in a condition ip* min < P + p*. If p* <C P, as one 
can expect for r*/ro ~ 1, we get for the boundary of the 
non growing phase 



p = 3p: 



(10) 



which, as Fig. [| demonstrates, is indeed in good agree- 
ment with the numerical data. 



C. Terrace formation. 

The average length T of a full swarming cycle was cal- 
culated by determining the position of the peak in the 
power spectrum of S(t) = J Q p*{x,t)dx, the time de- 
pendence of the total number of swarmer cells in the 
system. As Fig. || demonstrates, the dimensionless cycle 
time values are widely spread between values of 3 and 
12. However, T is only sensitive to changes in ro, r* and 
Pminl ' Pmin-, hence it does not depend on p max or r. 



5 




BB2000 



BB2235 t 9 



PRM1 



0.0001 0.001 , , 0.01 0.1 1 

V Pmin 

to Pmin 

FIG. 5. The cycle time T as a function of the approximate 
consolidation rate r* p* min . The data collapse indicates that 
the impact on T of the other parameters {p m ax and r) is neg- 
ligible. The various symbols correspond to different values of 
r*/r as 0.1 (+), 0.2 (x), 0.3 (□), 0.4 (o), 0.5 (A), 0.7 (v), 
1.0 (o). The dashed line is a plot of Eq.(|12[). Circles mark out 
the assumed parameter values characteristic of four different 
P. mirabilis strains. 

The average expansion speed v and terrace size w were 
also calculated in the parameter regime resulting oscilla- 
tory expansion of the colony. First we determined the 
time fi/3 when the system reached 1/3 of its maximal 
simulated expansion R m ax — R{t m ax), with R(t) be- 
ing the position of the expanding colony edge and t rnax 
is the total duration of the simulation. The average 
speed was then calculated for the time interval between 
Unin = max(*i/3> ~tmax - 5T) and tmax- for the last 5T 
long time interval, or for the last 2/3rd of the total expan- 
sion, depending on which was smaller. After obtaining 
v as [Rmax ~ R(t m in)]/(t max - t min ), the average ter- 
race width was calculated as w = vT. Fig. ||. shows the 
dependence of these parameters on the swarmer cell pro- 
duction density P and migration density threshold p* nin . 
In general, decreasing p* min or increasing P results in an 
increase in both w and v. As Fig. |?]. demonstrates, for a 
given r*, the relevant parameter controlling w is Pj p min . 

The results on the cycle time T (Fig. |[) can be in- 
terpreted as follows. As p* min is a good estimate on the 
density of swarmer cells in the expanding front, at the 
end of migration phase the vegetative cell density within 
the new terrace is given by r* p min T* with T* being the 
duration of the migration phase. Now the length of the 
consolidation phase, T~ T*, is determined by the require- 
ment that p must reach p m m- 



* o r (T-T*) 



o r T 



As T* 



X/tq, the estimate ( pT[ ) is simplified to 

i pmin*' 



r T = In ■ 



r* p 



(11) 



(12) 



which gives a rather accurate fit to the numerically de- 
termined data (Fig. ||). 




FIG. 6. Terrace size (a) and expansion speed (b) of 
the system for r*/r = 1 as a function of the dimension- 
less swarmer cell production density P. The connected 
points correspond to various values of the migration thresh- 
old Pmin/ Pmin- In general, decreasing p* min or increasing P 
results in an increase in both w and v. Only parameters re- 
sulting in periodic expansion were investigated. Circles mark 
out the assumed parameter values characteristic of four dif- 
ferent P. mirabilis strains. 
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FIG. 7. Terrace size w vs Pj p* min for various values of 
r*/ro. Note the collapse of the data presented in Fig. 6a. 
The terrace size vanishes approaching the parameter regime 
where no sustainable expansion of the system is possible. 
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FIG. 8. (a) Length Tl of the lag phase vs the initial inocu- 
lum density po for — 0, rotd = 3 and various values of r/ro 
(0.3, 0.1, 0.03, 0.01, 0.003 and 0.001). The transition between 
the regimes separated by p° m i n becomes more smooth for de- 
creasing r/ro. (b) The experimentally calculated tqTl values 
vs the estimated inoculum density (based on 5 mm and 7 mm 
inoculum droplet sizes for the RPM and BB strains, respec- 
tively) and the corresponding fits using Eqs. (|13[) and 



Since the duration (Tl) of the lag phase (the time pe- 
riod before the first migration phase) has been in the 
focus of many recent experiments, now we turn our at- 
tention towards this quantity. At least four processes 
determine Tl- First, there is a time ti associated with 
the biochemical changes required to switch into swarm- 
ing mode. As discussed in Sec. II. a, these processes 
take place only prior the first swarming phase, and are 
presumably related to sensing the altered environmen- 
tal conditions. Second, the cell population must reach 
the threshold density p° min (at time to). Third, td time 
is required to produce fully differentiated swarmer cells 
(at time ti), and finally, the density of the swarmer cells 
must reach the migration threshold p* m i n - Let us inves- 
tigate how these parameters depend on the initial inocu- 



lum density po- 

As p grows with a rate r$ until the appearance of 
swarmer cells, 



r h = max [t e + t d , h\(p mm /p )\ 



(13) 



The time development of p* can be estimated by the in- 
tegration of Eq.(|j|) with 7-* = (i.e., assuming rcf) 
yielding 



p*(t + tt) = / 9 e rotl [e (r °- r) ' - 1]. 

r -r 

Therefore, t* = Tl — t\ is given by 
1 



r t* 



1 - r/ro 



In 



ro-r p n 



Poe 



r h 



+ 1 



(14) 



(15) 



an expression usually giving a minor correction to t\. 

Fig. ||a. shows the above calculated toTl vs po for 
ti = 0. The increase in length of the swarmer cells was 
assumed to be 20- fold, thus rotd = In 20 « 3, which value 
can be seen for po S> Pmin- In the opposite limit, when 
Po *C Pmim we have tqTl ~ — lnpo + const. These rela- 
tions allow the determination of both r and r (using the 
known value of Pmin/Pmin) fr° m the experimental data 
on T L (po). 



E. Comparison with experiments 

Most of the published experimental data are related to 
the average period length, T and terrace size w. From 
these parameters the average expansion speed can be cal- 
culated as v = w/T, i.e., v is not an independent quan- 
tity. As we could see in the previous paragraph, from the 
density dependence of the lag phase the parameters ro, 
ti + td and r can be estimated. Notice that this estimate 
on ro is in principle different from the value obtained by 
the usual methods based on densitometry in liquid cul- 
tures. Technically, p* min could be also determined fL4[ , 
but such measurements are not published yet. 

There are four Proteus strains studied systematically 
in experiments: the PRM1, PRM2, BB2000 and BB2235 
strains (see Table I). To extract the values of the model's 
parameters the following procedure was applied, (i) We 
estimated r based on lag phase length measurements, 
(ii) From the calculated roT values the p* min l Pmin ratio 
was estimated (assuming r*/r — 1) based on Eq. (|l2|) , 
see Fig. [|. (iii) Using Eqs. (|l^ ) and (|l^), by a non- 
linear fitting procedure (Levenberg-Marquardt method, 
fL8| , see Fig. ||b.) p m in, r, and tt + td was determined. 
(The latter value is not relevant in respect the periodic- 
ity of the behavior.) (iv) Knowing p* min and p m im from 
the experimental terrace width data xq and P can be 
estimated using Fig. |. (v) Finally, t is given by 



tqt = In 



Pro 



(16) 
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The parameter values of the model are summarized in 
Table II, together with the predictions on T, v and w. 
An excellent agreement can be achieved with biologically 
relevant parameter values. 

Two classes of model parameters should be distin- 
guished: (a) the ones which are related to the growth 
and differentiation of the cells (tq, r* , p m i n , and P{r, r)) 
and (b) those which depend on agar softness (Dq and 
Pmin)- For a given strain we expect that a change in 
the agar concentration influences only the latter group, 
while changes in temperature may affect both, but pri- 
marily fQ. In fact, as Table I. demonstrates, by changing 
Dq and keeping all the other, growth-related parameters 
constant, we could quantitatively reproduce the colony 
behavior observed on various agar concentrations. Sim- 
ilar statement holds for the temperature effects as well, 
where the only parameter we changed was the growth 
rate tq. 



V. DISCUSSION 

Periodic bacterial growth patterns have been in the fo- 
cus of research in the last few years. Since a colony can 
be viewed as a system where diffusing nutrients are con- 
verted into diffusing bacteria, one may not be surprised 
by the emergence of spatial structures p9| . However, the 
periodic patterns of bacterial colonies are qualitatively 
different from the Liesegang rings (for a recent review 
see [M) developing in reaction-diffusion systems: the 
spacing between the densely populated areas is uniform 
and independent of the concentration of the other dif- 
fusing species, i.e., the nutrients. The Turing instability 
is also well-known for producing spatial structures pi) , 
but in that case the pattern emerges simultaneously in 
the whole system. It is also known, that bacteria can 
aggregate in steady concentric ring structures as a con- 
sequence of chemotactic interactions 22 2^|, but as we 
discussed in Sec. II., it is established that swarming of 
P. Mirabilis does not involve chemotaxis communication. 
Therefore, none of the well known generic pattern form- 
ing schemes can explain the colony structure of swarming 
bacteria. 

As we mentioned in the introduction, oscillatory 
growth is also exhibited by other bacterial species. One 
of them, Bacillus subtilis, has been the subject of system- 
atic studies on colony formation and a number of models 
have been constructed to explain the observed morphol- 
ogy diagram (for recent reviews see 
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. Only one 

model addressed the problem of migration and consolida- 
tion phases: Mimura et al |2q] set up a reaction-diffusion 
system in which the decay rate of the bacteria was depen- 
dent both on their concentration and the locally available 
amount of nutrients. The periodic behavior is then a con- 
sequence of the following cycle: if nutrients are used up 



locally, then the bacterial density starts to decay pre- 
venting the further expansion of the colony. Nutrients 
diffuse to the colony and accumulate due to the reduced 
consumption of the already decreased population. The 
increased nutrient concentration gradually allows the in- 
crease in population density and the expansion of the 
colony, which starts the cycle from the beginning. While 
this can be a sound explanation for B. subtilis, as we dis- 
cussed in Sec. II., the nutrient limitation clearly can not 
explain neither the differentiation nor the consolidation 
of P. Mirabilis swarmer cells. 

Another recent study |27j] focused on the swarming of 
Serratia liquefaciens. In that case the structure of the 
molecular feedback loops are better explored, and were 
resolved in the model. The production of a wetting agent 
was initiated by high concentrations of specific signalling 
molecules. The colony expansion was considered to be a 
direct consequence of the flow of the wetting fluid film, 
in which process the only effect of bacteria (besides the 
aforementioned production) was changing the effective 
viscosity of the fluid. The wetting agent production was 
downregulated through a negative feedback loop involv- 
ing swarmer cell differentiation. This scenario is certainly 
not applicable to P. Mirabilis, where swarmer cells ac- 
tively migrate outwards and their role is quite the oppo- 
site: enhancing the expansion of the colony. 

The first theoretical analysis focusing on P. mirabilis 
was performed by Esipov and Shapiro (ES) in p8[ . Their 
model was constructed based on assumptions similar to 
ours, and could reproduce the alternating migration and 
consolidation phases during the colony expansion. How- 
ever, the complexity of the ES model involves a rather 
large number of model parameters, which practically im- 
pedes both the full mapping of the parameter space and 
the quantitative comparison of the model results with ex- 
perimental findings. The major differences between our 
and the ES model can be summarized as follows: (i) we 
do not resolve the age of the swarmer population. In- 
stead, we have a density measure and a constant decay 
rate implying an exponential lifetime distribution on the 
(unresolved) level of individual cells. Since the available 
microbiological observations [p|Jf3| suggest only that the 
lifetime is finite, there is no reason for preferring any 
specific distribution, (ii) We did not incorporate into 
our model an unspecified "memory field' with a built-in 
hysteresis. Instead, we implemented a density- dependent 
motility of the swarmer cells, which behavior has been 
indeed observed . (iii) In our model the fully differ- 
entiated swarmer cells do not grow, which assumption is 
probably not fundamental for the reported behavior, but 
it seems to be more realistic because of the repression of 
many biosynthetic pathways ||]. Finally, (iv) we do not 
consider any specific interaction between the motility of 
swarmer cells and the non-motile vegetative cell popula- 
tion. Although such interactions probably exist, they are 
undocumented, and as we demonstrated, are not required 
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for the formation of periodic swarming cycles. However, 
such effects can be important in the actual determination 
of the density profiles. 

With these differences, which are not compromising 
the biological relevance of the model, we were able to 
map completely the phase diagram, establish approxi- 
mate analytical formulas and estimate the value of all 
model parameters in the case of four different strains. 
In addition, experimental data measured under various 
conditions could be explained with one particular param- 
eter setting in the case of the PRM1 strain indicating 
the predictive power of our approach. Our model is a 
minimal model in the sense that all of the explicitly con- 
sidered effects (thresholds, diffusion, etc.) were required 
to produce the oscillatory behavior, thus, it can not be 
simplified further. Such minimal models can serve as 
a comparison baseline for later investigations of various 
specific interactions. 

The values of the microscopic parameters of the model 
can be either measured directly (like r , P* m im Pmin, Pmax 
or r*) or can be determined indirectly from experimen- 
tal data (as p° min and r). Most of these measurements 
have not yet been performed, we hope that our work will 
motivate such experiments further examining the valid- 
ity of our assumptions. In fact, one of the parameters, 
r*/ro was set to 1 during the fitting processes, as cur- 
rently there is no available data to estimate its value. 
Our numerical results suggest that it is probably larger 
than 0.3, and it is unlikely to be larger than 2 (meaning 
an average lifetime less than 30 minutes). Within this 
range our qualitative conclusions are valid, while the nu- 
merical values of the parameter estimates can change up 
to a factor of 3. 

The behavior of "precocious" swarming mutants re- 
ported in H deserves special attention. First, we would 
like to comment on the huge difference found in the value 
of the transition rate r (see Table II). We emphasize that 
this is not an arbitrary output of a multiparameter fit- 
ting process. First, we have reasons to believe, that the 
motility thresholds of the two BB strains are rather sim- 
ilar. Knowing the growth rates and the cycle times, Eq. 
( |l2| ) shows us that the difference in the values of p m in 
can not exceed one order of magnitude. Assuming then 
this maximal difference in p m im T remains the only free 
variable in Eqs. (13)-(^5[), and the fitting can be per- 
formed unambigously. Thus, we are quite confident that 
such a large difference exists in r showing that the rsbA 
gene (in which these strains differ) influences not only the 
cell density threshold, but the rate of differentiation as 
well. It is also interesting to note that in Fig. 8b the be- 
havior of the PRM2 and PRM1 strains reflect a relation 
very similar to that of the BB2000 and BB2235 strains. 
Finally, our calculations predicted a slightly longer cycle 
time for the precocious swarming mutant BB2235, which 



is also in accord with the actual experimental findings 
(see Fig 2. of §). 

In our model the assumed functional form of the 
density-dependence of the diffusion coefficient is some- 
what different from the most often considered one p4j,p6l , 
namely 



D(p) ~ p k 



(17) 



The advantage of (O) is that it allows analytic solutions 
for certain cases pjE however, it describes an unlimited, 
arbitrarily fast diffusion inside the colony where the den- 
sity is high. In contrast, in real colonics the diffusion 
of cells is certainly bounded, and the expansion of the 
boundary can be often limited by the supply of cells from 
behind [Q. Therefore we believe that our thresholdcd 
formulation (Q) is a better approximation of what is tak- 
ing place inside the real colonies. 

Finally we would like to comment on the role of nu- 
triens in the swarming behavior of P. mirabilis. In 
our model there is a phenomenological parameter (r) 
determining how long the swarmer cells are produced 
at a given position in the colony. When investigating 
the dependence of the cycle time on this parameter, as 
Fig. 5 demonstrates, we found an extremely weak ef- 
fect. Thus, at least within the framework of this model 
there is no contradiction between the assumption that the 
swarmer cell production ceases due to nutrient (or accu- 
mulated waste) limitations, and the seemingly nutrient- 
independent cyclic behavior. In fact, this idea can be 
developed further. By increasing r (or decreasing the 
motility threshold p* min ) we arrive into a regime where 
the migration/consolidation phases are not clearly sep- 
arable as a motile swarmer cell population exists even 
when the expansion of the colony is slower. Experiments 
mapping the morphology diagram of P. mirabilis (Fig. 2 
of 1 29 1) showed that there are certain values of agar hard- 
ness and nutrient concentration, for which the expansion 
of the colony is still oscillating, but the periodic density 
changes are smeared out due to the presence of motile 
swarmer cells in the consolidation periods. If one asso- 
ciates the increasing agar hardness with increasing p* min 
and the nutrient concentration with r then one can qual- 
itatively reproduce those (i.e., the P r and Ph) regions of 
the morphology diagram. 



ACKNOWLEDGMENTS 

One of the authors (M.M.) is grateful to T. Matsuyama 
and H. Itoh for many stimulating discussions on experi- 
mental results. This work was supported by funds OTKA 
T0I9299, F026645; FKFP 0203/197 and by grants No 
09640471 and 11214205 from the Ministry of Education, 
Science and Culture of Japan. 



9 



Strain 




PRM1 


r>TJ TV /TO 


rJrJzUUU 




Experimental 


Temperature 




oZ O 










07OP 


07OP 


condition 


Agar 


2.0% 


2.45% 


2.0% 


n.a. 


n.a. 


2.0% 


n.a. 


n.a. 




Reference 














R 


R 


i 


1 




T[h] 


4.7 


4.7 


4.0 


3.5 


8.5 


6.0 


3.0 


3.1 


Colony-level 


v [mm/h] 


1.7 


0.6 


1.0 


n.a. 


n.a. 


n.a. 


3.3 


3.3 




w [mm] 


8.0 


3.0 


3.8 


n.a. 


n.a. 


n.a. 


10 


10 


Cellular-level 


TO [1/h] 


0.6 


0.6 


n.a. 


1.0 


0.4 


n.a. 


n.a. 


n.a. 



TABLE I. Summary of experimental data for four different P. mirabilis strains under various experimental conditions. The 
value of ro was determined from growth monitoring in liquid cultures. 



Strain 




PRM1 


PRM2 


BB2000 


BB2235 


Experimental 


Temperature 


32° C 


32°C 


32°C 


37°C 


22°C 


32°C 


37° C 


37°C 


condition 


Agar 


2.0% 


2.45% 


2.0% 


n.a. 


n.a. 


2.0% 


n.a. 


n.a. 


Microscopic 


ro [1/h] 


0.53* (0.6) 


0.53* (0.6) 


0.7* 


1.0 


0.4 


1.0* 


2.5* 


1.5* 


parameters 


p min [cell//um 2 ] 






0.06 






0.6 


2.0 


0.2 


(independent) 


Pmin [cell/Vm 2 ] 






6 ■ 10~ 3 






6 ■ 10" 3 


4 • 10~ 3 


4- 10" 3 




Pmax [cell//um 2 ] 






0.6 






240 


3000 


2.2 




r/r 






10- 1 






10" 4 


2 ■ 10~ 5 


3- 10" 2 




D [mm 2 /h] 


20 


3.2 


6 








60 


40 


(derived) 


Pmm [cell/Vm 2 ] 






3 • 10~ 3 






3 • 10" 2 


HT 1 


10" 2 




P 1 pmin 






0.9 






0.04 


0.03 


0.3 




. T/T 






0.7 






0.9 


0.9 


0.3 




Pmin/ Pmin 






10- 1 






10" 2 


2 • 10" 3 


2 • 10" 2 




xo [mm] 


6 


2.3 


3.0 








5.0 


5.0 




Wo [mm/h] 


50 


20 


27 








85 


70 


Macroscopic 


T[h] 


5.5 (4.7) 


5.5 (4.7) 


4.7 (4.0) 


3.3 (3.5) 


8.2 (8.5) 


5.6 (6.0) 


2.9 (3.0) 


3.2 (3.1) 


behavior 


v [mm/h] 


1.4 (1.7) 


0.5 (0.6) 


0.8 (1.0) 


n.a. 


n.a. 


n.a. 


3.4 (3.3) 


3.1 (3.3) 




w [mm] 


7.8 (8.0) 


3.0 (3.0) 


3.9 (3.8) 


n.a. 


n.a. 


n.a. 


10 (10) 


10 (10) 



TABLE II. Model parameters and the corresponding results for the strains and experimental conditions specified in Table 
I. The model has seven microscopic parameters, the rates ro, r, r* , the threshold densities pmin, pmin an( l pmax and the 
diffusivity Do- For each of the strains r* = ro was assumed. For comparison, other (derived) microscopic parameters are also 
included. The calculated period lengths, terrace sizes and expansion speeds are also presented together with the corresponding 
experimental values (in parentheses). The values marked by an asterisk (*) were derived from the lag phase length data based 
on Eqs. (|l^) and (jl|). Note the similarity between the PRM1 and BB2235, and also between the PRM2 and BB2000 strains. 
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